N = [2:65,70:5:200]';
%N = [2:65]';

x = load('tableData.dat');

Hnumerical = x(:,2) ./ 2;

o = optimset('MaxFunEvals',4000,'MaxIter',4000);

x1 = lsqcurvefit(@fitFunctionForH,[0,0,0,0,0,0],N,Hnumerical,[],[],o);
x2 = lsqcurvefit(@fitFunctionForHApprox,[0,0,0,0,0,0],N,Hnumerical,[],[],o);


plot(N,abs(fitFunctionForH(x1,N) - Hnumerical),'*',N,abs(fitFunctionForHApprox(x2,N)-Hnumerical),'.');